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Алгоритм решения краевых задач в системах 
компьютерной алгебры по *-методу Ланцоша 


В статье построен алгоритм т-метода Ланцоша для решения многоточечных линейных краевых задач 
для линейных дифференциальных уравнений порядка К с многочленными коэффициентами. По этому 
алгоритму в компьютерных системах символьного преобразования вычисляют многочлен порядка и. 
Доказана эквивалентность этого алгоритма и алгоритма применения а-метода Дзядыка. Результаты 
исследования а-метода доказывают существование решения исходной задачи по алгоритму, сходимость 
последовательности таких решений (с ростом параметра и алгоритма) к точному решению краевой 
задачи и точные и конструктивные априорные и апостериорные оценки погрешности в пространствах 
Сеь Е. [‹.5] ДлЯ достаточно широкого класса уравнений и краевых условий. 


Введение 


Задача. Построить алгоритм решения многоточечной линейной краевой задачи 
для ЛДУМК (®@$К) на отрезке [а, 6] в системах компьютерной алгебры (СКА). По этому 
алгоритму вычисляют решение, оптимальное для символьных преобразований в СКА. 

Многоточечная линейная краевая задача для ЛДУМК - это система, состоящая 
из линейного дифференциального уравнения с многочленными коэффициентами [1] 


БЫ + С= 0, РУ =Ау®+... + Су (1) 
и А линейных краевых условий в точках 4 1 = 4,1 = Г...., К, отрезка [а, Б] 
Сопа(у) = {П1у] = в; п, а 1=а; Е [а, 6], 1 =1,.... К}; (2) 


где О]у| - линейный дифференциальный оператор с постоянными коэффициентами 
и 
БУ] = А; У +...+ С;у 


имеет порядок К 1 = К; = о’ еди(Б1у]) < К = о еди(О)). Точки задания условий (2) 
являются не особыми точками ЛДУМК (1) - А(4;) #0, 1 = [,..., К. 

Актуальность задачи. Маре, Мафетайса, Ма@саа, МаЧаь, АР$ и другие компью- 
терные системы стали естественной средой математического моделирования. Эти систе- 
мы символьно преобразуют композиции специальных математических функций (СМФ) и 
решают дифференциальные уравнения — являются СКА. Краевые задачи (1), (2) являются 
классическим аппаратом математического моделирования [1]. 

Возможности СКА. СКА не имеют процедур для решения отдельных классов 
краевых задач. Система Маре вычисляет аналитическое решение обыкновенных диф- 
ференциальных уравнений в следующем виде: 

— композиция СМФ (она существует в очень редких случаях), 

— частная сумма ряда Тейлора решения задачи Коши порядка и-. Порядок этого много- 
члена, как правило, меньше 10. Ряд Тейлора является аппаратом аппроксимации функций, 
не оптимальным по точности. Следовательно, частная сумма ряда Тейлора решения 
задачи Коши для ЛДУМК (1) (вычисленная системой Мар!) часто не удовлетворяет 
основное требование математического моделирования — точность. 
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Алгоритм [2] решает класс краевых задач для ЛДУМК (1), (2). Для краевых условий 
этого класса задач существует многочлен порядка А - 1 


Ув! = 50е(Соп4(ук1 Е Рь 1) (3) 
и оператор 
К[и] = зое(у® = и, {РДУ] ха. =0,1=1,...,}), Рь — Ро (4) 
линейного преобразования функции и в первообразную этой функции порядка А 
Ки] укги =У® >, Свы > Сы. 
Эта первообразная удовлетворяет краевые условия (2). Многочлен ур! Е Рё! 
уравнений Соп(у,, = Рь)) имеет порядок К — 1, символьные коэффициенты и является 


элементом общего вида множества алгебраических многочленов порядка К - [. 
По алгоритму [2] вычисляют алгебраический многочлен порядка и 


Уи= с +сх+ ...+ сих" = авогийт (а а5К (1), (2), [а,6], п). (5) 
Коэффициент оптимальности алгоритма [2] в пространстве САаы ограничен 
Ск(айвотийт , паз4 (1), (2), Счаы ) = у — У» Почдаы / Е»С», С\ины) = О) (6) 
В коэффициенте оптимальности (6) эталоном для нормы погрешности алгоритма на 
краевой задаче (1), (2)) в пространстве Саы = С^К [а,Ь] 
|[у— у» см ьы = тах {|| у — У» [св ь | 7 — Ув) Пакьы „= | 0 Уна} 


является величина наилучшего приближения точного решения у исходной задачи 
1а5К алгебраическими многочленами порядка п в этом пространстве 


К р 
Ещу, СЧа.ы) = ШЁс 0, ... еп} [| У(®) - (со + ... + с,х) сан. 


Недостаток алгоритма [2]. СКА не преобразуют условия (2) в оператор Ки] (4). 
Метод решения задачи. //остроить алгоритм т-метода Ланцоша [1], эквивалент- 
ный алгоритму [2]. 


| Алгоритм 1 


Вход: О[у] + С =0 (1), Сопа(») (2), [а,Б]|, п. 
Выход: Многочлен (5). 


Преобразования: 
1. Вычислить аппроксимацию ЛДУМК (1) 
Бу, е Р„| + С + (Е и(2(х)) Е Ни) = 0. (7) 
2. Вычислить аппроксимацию краевых условий (2) 
Сопа(у, е Р„) = {1 у, Е Р |=: + :=0,1=1,...}. (8) 


3. Вычислить решение системы уравнений (7), (8) — многочлены 
{у, Етр(2(х))} = зо№е(О[у, Е Р„|+С + (Е„(2(х)) Е Ни») = 0, Сопа(у е Р„). (9) 


В уравнениях (7) и (8) многочлен у, е Р, имеет вид (5) и символьные коэффициенты. 
Поэтому он является элементом общего вида множества Р, — многочленов порядка и. 
Многочлен Ё„,„(2(х)) е Ни» уравнения (7) является дополнительным многочленом 
а-метода Дзядыка [3] с параметром р решения линейных интегральных уравнений 
Вольтерра с многочленными коэффициентами. Этот многочлен имеет параметры: 


т = 4е(Б[ у, ЕР,]+ 0), (10) 
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р=п-№, (11) 
где А = о’Я еди(Б[у]) - порядок старшей производной оператора О[}]|, вид 
(Еир(2(х)) е Нитр) — 5иБ5(х № 2(%), Етр), (12) 
где 
2(х) =2 (х-а)/(6Б-а-1, (13) 
Етр = 11 сйеБ(р + 1,х) + ... + тир сйеБ(т, х) (14) 


— дополнительный многочлен а-метода Дзядыка на отрезке [-1, 1]. Базис этого многочлена 
является ортогональным базисом пространства Гильберта [.2(а,Б; р) 


среь(О, 2(х)), сйеБ( 1, 2(х)), ..., спеБ(Е, х) = соз(1 агссо$(хХ)). (15) 


Многочлен Е„р(2(х)) е Нир имеет символьные коэффициенты (7,,..., р) = (ТКп),..., 
Ти-(п)). Поэтому он является многочленом общего вида пространства 
Нп\р — линейной оболочки элементов базиса (15) с индексом 1 =р + 1,..., т. 


2 Алгоритм 1 и классические методы 


Модельная краевая задача метода сплайнов [4, с. 63] 


у" = 100у, у(0) = 1, (О =1 (16) 
имеет единственное решение 
У=К[ и] +1 = с6(10х- 5) / сВ(5), и=у" = 100у, (17) 


где 
Ки] =“ а-хти@ &+'а-дхи@) 4. 


Функция у (17) — целая. Отличные от нуля коэффициенты Фурье — Чебышева функции 
у (17) на отрезке [0,1] только четные 


{аз Ку, [0, ])} 1 то = {0.47, 0.14, 0.021, 2 Я 
5.4 10°, 1.7 1077, 4.3 107 ‚8.5 101; 1.4 1012, 1.8 10-4}. (18) 


С ростом параметра { они регулярно убывают | а› у , [0,1]) | = 0(4)), 4 = 0.1. Следова- 
тельно, справедливо тождество 
Е у, Сю ) = (1 + 0(1)) | 42 изу у, [0, 1]) |. (19) 
Вычислительный эксперимент с алгоритмом 1. Решение краевой задачи (16) 
по алгоритму 1 на отрезке [0,1] является четным многочленом и справедливо тождество 
у» = а[вогийт_1(1а5К (16), [0, |, п) = Ури, 


Этот многочлен тождественен решению задачи (16) по алгоритму [2]. Норма 
погрешности решения краевой задачи (16) по алгоритму 1 в пространстве Со] 
убывает с ростом параметра п алгоритма 


{Пу- уз [сю |У- уз [ао иг’ = {0.24, 0.13, 0.011, 
0.0008, 2.8 10, 8.8 10’, 2.1 10%, 4.1 10°, 6.4 1023. (20) 


Из тождеств (18) -— (20) можно сделать следующий вывод. 
Вывод 1. Коэффициент оптимальности решения по алгоритму [ краевой за- 
дачи (16) в пространстве Сул] ограничен 
Сиавотийт_ 1, а5К (16), Сю, ) = (1+ о(1)) (5 + 922], (21) 
о? = 03, 1.3, 0.9, 0.4, 1.1, 0, —0.1, —0.2, —0.4}. 
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Норма погрешности решения краевой задачи (16) по алгоритму 1 в прост- 
ранстве С’, !] Убывает с ростом параметра и алгоритма 


[У Узнз [с ол = | Узна |е^холу = | У" — Уна ао = || — из [соль 
{Пи - м2; [со1}:=0 = 186, 22, 2.8, 0.24, 1.4 107, 
5.8 10", 1.8 10°, 4.5 107’, 8.8 109}. (22) 


Отличные от нуля коэффициенты Фурье — Чебышева прозводных у' и у" =и= 100 у 
функции у (17) на отрезке [0,1], с ростом параметра и, регулярно убывают (18). 
Следовательно, имеет место цепочка тождеств: 


Ежу, Сю) = Е»-у”, Сю) = (1 + 0(1)) [ аку", [0,1]) | = (100 + о(1)) | >», [0,1] | 


и из тождеств (18), (22) можно сделать следующий вывод. 
Вывод 2. Коэффициент оптимальности решения по алгоритму 1 краевой за- 
дачи (16) в пространстве Сул монотонно убывает к единице (с ростом параметра п ) 


С‚(авогийт_1, пак (16), Сол) = (1+ (1) (1 + а 2}, (23) 
{02-1 = {0.8, 0.6, 0.3, 0.2, 0.1, 0.08, 0.06, 0.05, 0.04, 0.03}. 


Коэффициент оптимальности (6) сплайн-метода [4]. По методу [4, с. 63] 
вычисляют сеточную функцию 


5 = {5-9 = зрИпе_тефоа(иазд (16), [0,1], 9) = но. 


Рассмотрим метод решения краевой задачи (16) — интерполяция Р. многочленом 
порядка 4 сеточной функции 5. (вычисленной сплайн-методом [4]) 


Ра5а] = РзрПпе теоа(а5к (16), [0,1], 4)]. (24) 


Норма погрешности этого метода в пространстве Сиы, очевидно, не меньше 
сеточной нормы погрешности сплайн-метода [4] 


плах, с [0,1] | 769) — Ра 54 || > тах, 41 УС ) — 34 И ) |. 


Поэтому значение коэффициента оптимальности метода Р.[54] (24) на краевой зада- 
че (16) не меньше 


С (Раба, (16), Сю, п) > тахго....а| У(ИЧ) — За(Ма) | / Е, Сю). 


Сравнение метода [4] и алгоритма 1. Погрешность аппроксимации точного 
решения (17) краевой задачи (16) сеточной функцией 5, принимает [4, с. 63] следующие 
значения: 


на 


Поэтому из тождеств (18), (19) можно сделать следующий вывод. 
Вывод 3. Коэффициент оптимальности метода (24) на задаче (16) растет 
экспотенциально (с ростом числа узлов сетки) 


Сто(Рю[5 16], (16), Сю, 1) > 2000, С›о(Р2 [$20], (16), Ср, 1) > 10'2/6, 
С(Ра 8], (16), Сю,н) = О(2"), => 2. 


Такая же зависимость коэффициента оптимальности метода от параметра 4 (порядка 
вычисляемого многочлена или числа узлов сетки) имеет место и в случае решения 
краевой задачи (16) по другим методам с насыщением. На основании этой зависимости и 
тождеств (21), (23) можно сделать следующее заключение. 

Вывод 4. Решение краевой задачи (16) по алгоритму [1 оптимально для символьных 
преобразований в СКА. (Это основное отличие алгоритма 1 от сплайн-метода [4, с. 63], 
метода (24) и других методов с насыщением.) 


«Штучний 1нтелект» 12008 41 


1д. Денисенко П.Н. 


3 Эквивалентность алгоритма 1 и алгоритма [2] 


Алгоритм [2] 
1. Вычислить интегральный оператор К[и] (4) и многочлен уд. (3). 
2. Преобразовать систему уравнений (1), (2) в линейное функциональное уравнение 


Ни] +8=0, Ми] = Б[КТи, & = Буый+ 6, и=у”. (25) 
3. Вычислить аппроксимацию уравнения (25) по а-методу [3] с параметром р=п-^ 
Мир Е Р» | +& + (Епр(20)) Е Ни) = 0, (26) 
где дополнительный многочлен имеет параметр 
т = де>([ и, Е Р}| + 8). (27) 


4. Решить уравнение (26) и вычислить многочлены и,, Еш,р(2(Х)). 
5. Вычислить преобразование К[и]+у- (3), (4) многочлена и, 


Уи = К[Ир] + У (28) 


— искомую аппроксимацию решения у краевой задачи (1), (2). 

Замечание 1. Согласно определению оператора К (4) и многочлена уз! (3), 
многочлен у, (28) является решением СЛАУ (8). 

Теорема 1. Пусть: . 
— оператор Ки] (4) преобразует моном х', Г Е М в полином 


КЕ Вх... +е, Ь,=0, (29) 
— параметр алгоритма 1 (и алгоритма [2]) 
п>2к, К= ога еда(О[у]|). (30) 


Тогда уравнение (26) тождественно замене в уравнении (7) многочлена уи Е Ри 
на многочлен Ки, ЕР, | + ул, где р =п-- уравнению 


$165((у, ЕР,) = Ки, ЕР» +уы, Б[у, ЕР] + С + (Ем (2х) Е Нть)=0) 
= (ЖК и, ЕР] + ук] + С + (Еир(2(х)) Е Нть), (31) 


где т = 4е(ОБ[у, ЕР,| +0). 

Доказательство. Согласно определению уравнения (25), первое слагаемое уравне- 
ния (7) тождественно первому слагаемому уравнения (26). 

Параметр р дополнительного многочлена ЁЕ»,,(2(х)) уравнений (7), (31) и (26) 
тождественен по определению. 

Тождество первого параметра т (10) и (27) дополнительных многочленов этих 
уравнений 

д4е2(Б[у, ЕР] + С) = 4е(Ци, ЕР, ] + 8), р=п-К 


в случае 4е2(Б[ у, е Р‚] + С) = 4е5(С) очевидно. 

Тождество первого параметра т (10) и (27) дополнительных многочленов этих 
уравнений в случае 4е(О[ у, е Р),]| + О) = 4е>(Б[]у, е Р)„|) непосредственно следует из 
очевидного тождества 


Чез(О[К[и» = Р»] + уь1]) = ае(О[К[и» = Р»]]), 
неравенств 
де(В[Ки, ЕР)}]]) < 4е>(Б[ у, ЕР„]), п =р- А, (32) 
де (В[Ки, е Р„]]) > тах {4ез(Б[ а д4е(В[х” ]}}, (33) 
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тождества 
4ез(О[», е Р„]) = тах{аез(РП]), 4ез(О[х]), ..., Че (Ох ])} (34) 
и тождества, имеющего место при выполнении условия теоремы и > 2, 


де®(О[у, е Р‚|) = тах {Чео(Р[х”), .... Че2(О[х"]}}. (35) 


Тождества (34), (35) доказаны нами в работе [5]. 
Лемма 1. //Лусть Б[Уу| — линейный дифференциальный оператор с многочленными 
коэффициентами вида (Т) и для оператора К|и] (4) справедливо тождество (29). 
Тогда справедливы неравенства (32), (33). 
Доказательство. Согласно определению (4) оператора Ки], многочлен Ки, Е Р}] 
имеет порядок 
4е2(К и» ЕР›„]) =р+к=и. 


Поэтому он является частным случаем многочлена у, Е Р,. Следовательно, много- 
член О[Ки, е Р)„]] является частным случаем многочлена О|[у, е Р‚,]| и порядок много- 
члена О[К]и, е Р),]] не превышает порядок многочлена Б[у, е Р,| — имеет место 
неравенство (32). 

Согласно тождеству (29), имеет место тождество 


Кс] = со К] = со (0х +...+ ео), Б0=0 
и неравенство (33) в случае р = 0 очевидно 


дез(О[К[со]]) = 4е®(Б[К1]) = тах {Че(Б[х“), ...} > 4е(Б^). 
Операторы О[у] и КГи] — линейные. Поэтому имеет место тождество 
БК и, е Р = Б[К[и,1 Е Р.И + © КИ. 


Согласно этому тождеству и правилам вычисления порядка многочлена, имеет 
место тождество 


де›(О[К]и; е Р;]]) = тах{4е(Б[Ки,-1 е Р..1 ||, Че( ВК. (36) 
Согласно тождеству (29), имеет место неравенство 
Чев(Р(КТх) > 4е(х"}. 

Если неравенство (33) имеет место дляр =5- 1 ‚то из этого неравенства и тождества (36) 
следует справедливость неравенства (33) в случаер =5. 


4 Сходимость алгоритма 1 и оценки погрешности 


Согласно теореме 1, для алгоритма 1 имеют место результаты [2] исследования 
алгоритма [2]. Следовательно, для достаточно широкого класса ЛДУМК (1) и краевых 
условий (2), эти результаты: 

— доказывают существование решения у/ (9) системы аппроксимирующих уравнений (7), 
(8) алгоритма 1; 

— доказывают сходимость последовательности многочленов уд ‚ Ук+1 ,... (9) к точному 
решению у краевой задачи для ЛДУМК (1), (2); 

— устанавливают точную и конструктивную апостериорную оценку нормы погрешности 
алгоритма 1 в пространстве Си,5] и Сы 

— устанавливают точную и конструктивную априорную оценку коэффициента оптималь- 
ности алгоритма 1 в пространстве С |и.5]; 

— доказывают ограниченность коэффициента оптимальности алгоритма 1 в прост- 
ранстве Саь. 
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5 Вычисление оценок оптимальности алгоритма 1 


Согласно теореме 1, алгоритмы [2] вычисляют: 
— точную оценку нормы погрешности алгоритма 1 в пространстве Ср] и Саь } 
— точную оценку коэффициента оптимальности алгоритма 1 в пространстве С’«.ьдля 
достаточно широкого класса ЛДУМК (1) и краевых условий (2). 


6 Программирование алгоритма 1 в СКА 


Операторы СКА не преобразуют систему уравнений (1), (2) в систему уравне- 
ний (7), (8) и не решают систему уравнений (7), (8). Поэтому мы построим алгоритм 2 
преобразования краевой задачи для ЛДУМК (1), (2) в искомый многочлен у, и докажем 
эквивалентность аппроксимации по этому алгоритму краевой задачи для ЛДУМК (1), (2) 
системе уравнений (7), (8). 
Алгоритм 2 
Вход: О[у] + С =0 (1), Сопа(») (2), [а, Б], п. 
Выход: Многочлен у„ (9). 
Преобразования: 
1. Вычислить оператор О[у] + С — левую часть ЛДУМК (1). 
2. Вычислить К — порядок старшей производной оператора О|[у]. 
3. Вычислить многочлен у, Е Р,. 
4. Вычислить многочлен О|[у, ЕР) + Ц. 
5. Вычислить порядок т (10) многочлена О[у, ЕР) + Ц. 
6. Вычислить параметр р (11) дополнительного многочлена. 
7. Вычислить многочлен Е» (14) (с символьными коэффициентами). 
8. Вычислить многочлен 2(х) (13). 
9. Вычислить многочлен Еш,р(2(х)) Е Ни} (12). 
10. Вычислить аппроксимацию оператора О[у] + С (1) - многочлен 


ОБ] у, ЕР, + С 1 (Епр(2(%)) Е Нт,. (37) 
11. Вычислить систему линейных алгебраических уравнений 
5 = {# (Оу, ЕР, +С+ (Е „р (20)) Е Нть)) =0, 1=0,..., т}, (38) 


где {(Р) — коэффициент монома порядка 1 полинома Р. 

12. Вычислить объединение СЛАУ (38) и СЛАУ (8). Для! =1,..., А вычислить: 
12.1. Условие (2) с номером [. 
12.2. Линейное дифференциальное уравнение этого условия 


БУ] + С; 5 0. (39) 


12.3. Левую часть уравнения (39) — оператор ОДу] + С;. 
12.4. Многочлен Эду, ЕР) + С:. 
12.5. Точку задания краевого условия 1 (2). 
12.6. Левую часть уравнения 1 СЛАУ (8) — (Ру, ЕР) + С) [=а +. 
12.7. Уравнение # СЛАУ (8) В 
5; = ((Р4у, ЕР + С) [5=а:= 0). (40) 


12.8. Объединение уравнения (40) с объединением СЛАУ (38) и вычисленной ранее 
частью СЛАУ (8) 


сопс(5;= 0, {5 (38), 5 =0, ..., 5 = 0}). 
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13. Решить объединение СЛАУ (8) и СЛАУ (38) — вычислить значения коэффициентов 
Со ›...› Спь 1 ь...› Пр МНОГОЧЛеНОВ Ур , Ёт,р(2(Х)) (9): 


Сое{ = 5о[№е(5 (38), (8)) = {соеКу»), соеКЕт„ь)}. (41) 


14. Вычислить многочлен у, (5) с числовыми коэффициентами — преобразовать значения 
коэффициентов этого многочлена соеКу„) = {4, ...,е} (41) в сумму мономов: 


у» = зе" (соеКу,)) = а-...+ех". (42) 


Доказательство. Операторы Оу] условий (2) — линейные. Поэтому левая часть 
уравнения 7 системы (8) имеет вид 


(Гу, ЕР) + С). -а:=А, (у Е Р-Р, -4:+... + С. бл еЕ Ра: + 6, 


и это уравнение (относительно коэффициентов многочлена у, е Р‚)-— линейное. 
Функционал 1; — линейный и оператор О[у] — линейный. Поэтому левая часть урав- 
нения 1 = 0, 1,..., т системы (38) имеет вид 


КО, е Р„]) + КС) - ЩЕ»(2(х)) Е Нтр), 
где слагаемые являются линейными формами коэффициентов многочленов у, Е Ри 
Етр(2(Х)) Е Нтр: 


(Оу Е Р„|) = КОП] со +...+ К ЬОГ”]) си ‚ КС) = Соп5ь 
КЕт.ь(2(х)) Е Нир) = (свеб(р-+1, 2(х)) пи +...+ Ксбеб(т, 2(х)) сир. 


Многочлен (37) тождественен левой части уравнения (7) по определению. Этот много- 
член имеет порядок т (10). Следовательно, уравнение (7) эквивалентно СЛАУ (38) и 
многочлен у, (42) тождественен решению у„ системы уравнений (7), (8). 


7 Программирование алгоритма 2 в системе АРЗ 


Структура данных на входе 
1. ЛДУМК (1) имеет вид 


Я о и 


где у — атом, коэффициенты А ,..., С, С являются многочленами переменной (атома) 
х (термами) и имеют вид, естественный для математики. 
2. Краевые условия (2) являются списком вида 


2 =, д а О уно, 


3. Отрезок аппроксимации [а ‚ 6] определяет список (а, Б). 
4. Параметр и алгоритма является целым числом. 
Структура данных на выходе 
Многочлен у» (42) имеет числовые коэффициенты и вид, естественный для 
математики 
о о ет 


АРГАМ№М-процедура. Алгебраическая спецификация алгоритма 2. 
1её (тромМк, Бу = 0); в, ео 


К := ока еаа (ру); /* порядок Бу */ 

уп := па1п_ро1 (п); ум дамр м № 

Би зе апрт Е ааа Ву а У 62 

п := ед (сапр1Е (е1п_ро1 (Рп))); /* аезч(Р[у в]+6) */ 


р := пП- К; 
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Ем := Ер1 (п, ш-о+К, К); в мы в 
Ь := ага (106екуа1, 2); а := ага (106екуа1, 1); 
В := Б -а; 
я += сайр1{((276) * (5 + 5-1) ха +1); 7820) * 
Ве» бапр Ева = м ва ее а 
Вп: ==> вапр1Е (Бр Ем); “у ем */ 
5 := ро1 еча(Рп, ш); /* СЛАУ - аппр. ЛДУМК */ 
/х АППРОКСИМАЦИЯ условий (2) */ 
Во К 1 Е, 


Сова еее (сона». 
СопЯ --> ага (Сопа, 2); 


ро => ата (соло; 1) дж =са 

УМЕ 1 ==> аза{сова ль и Ге, = 0 

1её (ТРОМК_1, Бу = 0); /* Пу = р Е У] + 61%*/ 

Ве я <апршЕкВию ох ву, ре Е вв 
] 


.) 
В ыса-а воре прощании 6 
3>==> с0пс{5, `сору(5: 1) =.0) 
); 
ро1пЕ --> ага (Сопа, 1); д: а, 7 
ББУМК 1 == ава (Сова, 2) /* р К[ у] = 0 */ 
1е- (ТРОМКк_1, Ру = 0); /* ру = В к у] + ск */ 
Вузе авео ое КГ рее 97 
Вт ее вое (роте, вв ка 1жа к СЕ 
Бе» Оно (9» сор (= ОЕ И ® СЛАУ (38): {8) 97 
Хп := с; СоеЕЁ := зо1уе($); /* решение СЛАУ */ 
У ева, Сов) /* аппроксимация у */ 


Структура выхода операторов АРГАМ-процедуры 
Многочлен у, е Р, АРГАМ-процедура вычисляет в виде 


© (0) + с. хх ео Хх т, 
где х и коэффициенты с (0) ,... Сс(п) являются атомами. Оператор О[у| (Т) преобра- 


зует этот многочлен в многочлен О[у, = Р‚| порядка т (10). Оператор сапр1Е АРГАМ- 
процедуры приводит этот многочлен к сумме слагаемых вида 


с (1) *х ^ ] $, #=0....п,7=0,...,т, (43) 
где $ — операция умножения терма на число. 
Дополнительный многочлен ЕР, (14) АРГАМ-процедура вычисляет в виде 
с(п + 1) * стеб (р + 1,х) +...+ с(м + К) * себ (и, х), 
где коэффициенты с (п + 1),.... с(м + К) являются атомами, и многочлены Чебы- 
шева первого рода — стер (1,х) — вычисляет внутренний оператор солвера е"_5о[уе.ехе 
системы АР. Этот многочлен и его преобразование (12) имеют порядок т (10). Оператор 


сапр1Е АРГАМ-процедуры приводит многочлен Е»„»(2(х)) = Ни) (12) к сумме сла- 
гаемых вида 


а а (44) 
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Поэтому объединение СЛАУ (8) и СЛАУ (38) является списком уравнений вида 


СУ 5. КЕ Ее ОО) ча О 
и решение (41) объединения СЛАУ (8) и СЛАУ (38) является списком тождеств 


с (0) = В ,..., с(п) =р,..., см + К) =а. 


8 Исследование АРГАК-процедуры 


Теорема 2. Пусть коэффициенты многочленов и другие параметры исходных 
данных АРГАМ№-процедуры являются целыми или рациональными числами. 

Тогда АРГАМ-процедура вычисляет многочлен у (42) точно. 

Доказательство. АРГАМ-процедура имеет известные операторы. Эти операторы в 
системе АР$ преобразуют рациональные числа (коэффициенты многочленов) точно — не 
вносят в результат преобразований погрешности от выполнения арифметических 
операций компьютером с фиксированной длиной машинного числа. 

Теорема 3. АРГАМ-процедура имеет по параметру п полиномиальную сложность 

(т + 1) О(сапр! т) + О(и^), т = де(Б[у, е Р]+С) =п+О0(1), (45) 
где О( сапр/, т) — сложность преобразования оператором сапрИ многочлена Р: 
многочлен сапрКР) является суммой мономов вида (43) и (44). 

Доказательство. АРГАМ-процедура линейная. Следовательно, вычислительная 
сложность этой процедуры тождественна сумме вычислительной сложности ее операто- 
ров. Операторы АРГАМ-процедуры, за исключением оператора сапр, имеют по пара- 
метру и полиномиальную сложность О(и“), 5 < 3. Наибольшую вычислительную 
сложность имеет оператор ро! еда. Вычислительная сложность оператора ро! еди 
тождественна оценке (45). Оператор сапр преобразует многочлены О[у, Е Р)]| + С, 
Пу, ЕР) + Сь (12) и (37) в сумму мономов вида (43) и (44). 


9 Вычислительный эксперимент с АРГАМ-процедурой 


Описание краевой задачи (16) для АРГАМ-процедуры 
ргосез$ [1] := ( 


ТРОМК := (А1Е(у, 2) + (-100) *у=О0); 
Ссора := ((х=0, у+-1=0), (х=1, у+-1=0)); 
1пфехуа1 := (0, 1); че 
Результат преобразования процедурой задачи (16) и параметра и =2 
у 5 ад. К ЕЕ ях. 
Заключение 


Построенный в работе алгоритм т-метода Ланцоша оптимально решает в СКА 
достаточно широкий класс многоточечных линейных краевых задач для линейных 
дифференциальных уравнений порядка К с многочленными коэффициентами (1), (2). 


Дополнение 


1. Алгоритм 2’. Модификация алгоритма 2 — дополнительный многочлен 
СЛАУ (38) имеет ортогональный базис пространства Гильберта Н. 
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2. Модификация АРГАМ-процедуры. Оператор вычисления многочлена 
Чебышева заменен оператором вычисления многочлена базиса алгоритма 2' на отрезке 
[-1,1. Эта модификация АРГАМ-процедуры реализует алгоритм 2" в системе АР5. 

3. Для алгоритма 2' справедлив аналог теоремы 1. В этом аналоге уравнение (26) 
имеет дополнительный многочлен с базисом пространства Н. 
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П.М. Денисенко 

Алгоритм розв?язання крайових задач в системах комп'ютернот алгебри за т-методом Ланцоша 
У статт! побудовано алгоритм т-методу Ланцоша для розв’язання багатоточкових лиийних крайових 
задач для лиййних диференщальних р1внянь порядку К з багаточленними коефицентами. За ним в 
комп’ютерних системах символьного перетворення обчислюють багаточлен порядку п. Доведена 
екв1валентн!сть цього алгоритму та алгоритму застосування а-методу Дзядика. Результати теорй 
а-методу доводять 1снування розв’язання вих1дно! задач! за цим алгоритмом, з9жнасть посл1довност! 
таких розв’язань (з ростом параметру п алгоритму) до точного розв’язання крайово! задач, точн! 1 
конструктивн! апр1орн! та апостерорн! ощнки похибки в просторах Сы Сон для досить широкого 
класу крайових задач. 


Р.М№. Беш5епко 

Тве А!согйт о? Зоушо Фе Воипдагу-уаше РгоШетз$ ш е Сотшршег А!оебга Зузетз$ зто бе 
Гапс70$ с-тевоЯ 

Тве Гапс7о$ т-тейфо4 а]еогибт оЁ зоуше фе ши йрош{ Боипдагу-уаме ргоет юг Ппеаг Чегепйа1 
едчайоп$ оЁ огаег К \Л роупопта| соеЁЯслетв 1$ 4еуеоре4 ш фе агасе. ТЬе ро[упопиа| оЁ отаег п 1$ 
сотршеа Бу 15 а1еогИ т ш Фе сотршег а|еебга зуз{етл$. \Уе ргоуе4 Фе едуаепсе оЁ 1$ а1еогИт {0 
Фе У. К. РиуадуК а-тео4 а1согибт. Те гезеагсВ гези$ о Фе а-те#о4 ргоуе Фе зоаНоп ех15епсе Рог 
Фе ша] рго ет Бу Фе а!еогИбт, Бе сопуегоепсе оЁ Фезе зо[аНопз зедиепсе (уу фе шсгеазе оЁ Ше 
а]еогИйт рагатеег и) ю Ше ехас+ Боипдагу-уаше рго ет зоайоп, бе ехасё ап4 сопзиийуе езИтайез а 
рог! ап4 а роменог! т ве зрасез Сьы, С\аы] 


Статья поступила в редакцию 25.07.2007. 
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